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Abstract. A major factor affecting the readability of a graph drawing is its 
resolution. In the graph drawing literature, the resolution of a drawing is ei- 
ther measured based on the angles formed by consecutive edges incident to a 
common node (angular resolution) or by the angles formed at edge crossings 
(crossing resolution). In this paper, we evaluate both by introducing the notion 
of "total resolution", that is, the minimum of the angular and crossing resolu- 
tion. To the best of our knowledge, this is the first time where the problem of 
maximizing the total resolution of a drawing is studied. 

The main contribution of the paper consists of drawings of asymptotically op- 
timal total resolution for complete graphs (circular drawings) and for complete 
bipartite graphs (2-layered drawings). In addition, we present and experimen- 
tally evaluate a force-directed based algorithm that constructs drawings of 
large total resolution. 

Date: September 10, 2010. 

1 Introduction 

Graphs are widely used to depict relations between objects. There exist several criteria 
that have been used to judge the quality of a graph drawing |3I16) . From a human 
point of view, it is necessary to obtain drawings that are easy-to-read, i.e., they should 
nicely convey the structure of the objects and their relationships. From an algorithmic 
point of view, the quality of a drawing is usually evaluated by some objective function 
and the main task is to determine a drawing that minimizes or maximizes the specific 
objective function. Various such functions have been studied by the graph drawing 
community, among them, the number of crossings among pairs of edges, the number 
of edge bends, the maximum edge length, the total area occupied by the drawing and 
so on. 

Over the last few decades, much research effort has been devoted to the problem 
of reducing the number of crossings. This is reasonable, since it is commonly accepted 
that edge crossings may negatively affect the quality of a drawing. Towards this di- 
rection, there also exist eye-tracking experiments that confirm the negative impact of 
edge crossings on the human understanding of a graph drawing [20121122] . However, 
the computational complexity of the edge crossing minimization problem, which is 
A/'P-complete in general |11| , implies that the computation of high-quality drawings 
of dense graph is difficult to achieve. 

Apart from the edge crossings, another undesired property that may negatively 
influence the readability of a drawing is the presence of edges that are too close 
to each other, especially if these edges are adjacent. Thus, maximizing the angles 



among incident edges becomes an important aesthetic criterion, since there is some 
correlation between the involved angles and the visual distinctiveness of the edges. 

Motivated by the cognitive experiments by Huang et al. |14I15| that indicate that 
the negative impact of an edge crossing is eliminated in the case where the crossing 
angle is greater than 70 degrees, we study a new graph drawing scenario in which both 
angular and crossing resolution^ are taken into account in order to produce a straight- 
line drawing of a given graph. To the best of our knowledge, this is the first attempt, 
where both angular and crossing resolution are combined to produce drawings. We 
prove that the classes of complete and complete bipartite graphs admit drawings 
that asymptotically maximize the minimum of the angular and crossing resolution 
(Section [3]). We also present a more practical, force-directed based algorithm that 
constructs drawings of large angular and crossing resolution (Section @| . 

1.1 Previous Work 

Formann et al. [3] were the first to study the angular resolution of straight-line draw- 
ings. They proved that deciding whether a graph of maximum degree d admits a 
drawing of angular resolution (i.e., the obvious upper bound) is NT '-hard. They 
also proved that several types of graphs of maximum degree d have angular resolution 
0(2)- Malitz and Papakostas [T5] proved that any planar graph of maximum degree 
d. admits a planar straight-line drawing with angular resolution f2(Jj). Garg and 
Tamassia |12| showed a continuous tradeoff between the area and the angular resolu- 
tion of planar straight-line drawings. Gutwenger and Mutzel |13| gave a linear time 
and space algorithm that constructs a planar polyline grid drawing of a connected 
planar graph with n vertices and maximum degree d on a (2n — 5) x (|n — |) grid with 
at most 5n — 15 bends and minimum angle greater than |. Bodlaender and Tel [5] 
showed that planar graphs with angular resolution at least ^ are rectilinear. Recently, 
Lin and Yen [17] presented a force-directed method based on edge-edge repulsion that 
leads to drawings with high angular resolution. In their work, pairs of edges incident 
to a common node are modeled as charged springs, that repel each other. 

A graph is called right angle crossing (or RAC for short) graph if it admits a 
polyline drawing in which every pair of crossing edges intersects at right angle. Didimo 
et al. [6] showed that any straight-line RAC drawing with n nodes has at most An— 10 
edges. Angelini et al. [T] showed that there are acyclic planar digraphs not admitting 
straight-line upward RAC drawings and that the corresponding decision problem is 
MV-haiA. They also constructed digraphs whose straight-line upward RAC drawings 
require exponential area. Di Giacomo et al. [5] presented tradeoffs between the crossing 
resolution, the maximum number of bends per edges and the area. Dujmovic et al. 
[7] studied a Angle Crossing (or a AC for short) graphs, that generalize the RAC 
graphs. A graph is called aAC if it admits a polyline drawing in which the smallest 
angle formed by an edge crossing is at least a. For this class of graphs, they proved 
upper and lower bounds for the number of edges. 

Force-directed methods are commonly used for drawing graphs [8110] , In such a 
framework, a graph is treated as a physical system with forces acting on it. Then, 
a good configuration or drawing can be obtained from an equilibrium state of the 

The term angular resolution denotes the smallest angle formed by two adjacent edges 
incident to a common node, whereas the term crossing resolution refers to the smallest 
angle formed by a pair of crossing edges. 
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system. An overview of force-directed methods and their variations can be found in 
the graph drawing books |3I16) . 



2 Preliminaries and Notation 

Let G = (V,E) be an undirected graph. Given a drawing r(G) of G, we denote by 
p u = (x u ,y u ) the position of node u £ V on the plane. The unit length vector from 
p u to p v is denoted, by p u pl, where u,v £ V. The degree of node u £ V is denoted 
by d(u). Let also d(G) = max ue y d(u) be the degree of the graph. 

Given a pair of points q± , q^ & ^- 2 , with a slight abuse of notation, we denote by 
\\qi — q2 1 1 the Euclidean distance between qi and q^. We refer to the line segment 
defined by q\ and 52 as q\qi. 

Let ~ct and 7 be two vectors. The vector which bisects the angle between ~tit and 
~zf is ypg^ + n^ii . We denote by Bsc(7^, 7 f) the corresponding unit length vector. 

Given a vector p , we refer to the unit length vector which is perpendicular to p 
and precedes it in the clockwise direction, as Perp(T^). Some of our proofs use the 
following elementary geometric properties: 



tan (loi — L02) 



tanwi — tancj2 



(1) 



1 + tancji ■ tanw2 

tanw > uj (3) 



tan (w/2) 



cos LJ 



3 Drawings with Optimal Total Resolution for Complete and 
Complete Bipartite Graphs 

In this section, we define the total resolution of a drawing and we present drawings of 
asymptotically optimal total resolution for complete graphs (circular drawings) and 
complete bipartite graphs (2-layered drawings). 

Definition 1. The total resolution of a drawing is defined as the minimum of its 
angular and crossing resolution. 

We first consider the case of complete graphs. Let K n = (V, E) be a complete 
graph, where V = {uq, u\, . . . , u n — 1} arid E = V x V. Our aim is to construct a 
circular drawing of K n of maximum total resolution. Our approach is constructive and 
common when dealing with complete graphs. A similar one has been given by Formann 
et al. [H] for obtaining optimal drawings of complete graphs, in terms of angular 
resolution. Consider a circle C of radius r c > centered at (0, 0) and circumscribe 
a regular n-polygon Q on C. In our construction, the nodes of K n coincide with the 
vertices of Q. W.l.o.g., we further assume that u\, U2, ■ ■ ■ ,u n appear in this order in 
the counter-clockwise direction around (0,0), as illustrated in FigfTal 

Theorem 1. A complete graph K n admits a drawing of total resolution 0(— ). 
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Proof. We prove that the angular resolution of the presented drawing of K n is 
whereas its crossing resolution is — . First, observe that the arc of circle C that con- 
nects two consecutive nodes m and Ufj+iWodn is equal to for each i = 0, 1, . . . ,n — 
1. Therefore, the angular resolution of the drawing is — , as desired. Let now e; = 
(ui,Ui') and ej = (v,j,Uji) be two crossing edges. Without loss of generality, we as- 
sume that i < j < i' < j', as in Fig[Tal The crossing of ej and ej defines two angles 
<f> c and 4>' c such that cf> c + <f>' c = n. In Fig[Ia] </> c is exterior to the triangle formed by 
the crossing of e, and ej and the nodes Uj and u%' (refer to the dark-gray triangle of 
FigOU). Therefore: 4> c = (j' + Similarly, # = (i'_j)i + („-(j'-i))2. 
In the case, where j = (i + 1) mod n and j' = (i' + 1) mod n (i.e., the nodes Ui (v,i>, 
resp.) and Uj (Uj>, resp.) are consecutive), the angle <f> c receives its minimum value, 
which is equal to =^-. Similarly, we can prove that the minimum value of c/>' c is also 
— . This establishes that the crossing resolution is — . □ 

We now proceed to consider the class of complete bipartite graphs. Since an n- 
vertex complete bipartite graph is a subgraph of a n- vertex complete graph, the bound 
of the total resolution of a complete bipartite graph can be implied by the bound of 
the complete graph. However, if the nodes of the graph must have integer coordinates, 
i.e., we restrict ourselves on grid drawings, few results are known regarding the area 
needed of such a drawing. An upper bound of 0(n 3 ) area can be implied by [2]. 
This motivates us to separately study the class of complete bipartite graph, since 
we can drastically improve this bound. Note that the tradeoff between resolution 
and area has been studied by the graph drawing community, in the past. Malitz 
and Papakostas |18) showed there exist graphs that always require exponential area 
for straight-line embeddings maintaining good angular resolution. The claim remains 
true, if circular arc edges are used instead of straight-line (4]. More recently, Angelini 
et al. [I] constructively showed that there exists graphs whose straight-line upward 
RAC drawings require exponential area. 

Again, we follow a constructive approach. First, we consider a square 1Z = ABB A 
where its top and bottom sides coincide with H\ and £2, respectively (see Fig llbl) . 
Let H be the height (and width) of 1Z. According to our approach, the nodes of 
Vi {V2, resp.) reside along side rA (AB, resp.) of 1Z. In order to specify the exact 
positions of the nodes u\, u\, ■ ■ ■ , along side rA, we first construct a bundle of 
m semi-lines, say ly,. . .,£ m , each of which emanates from vertex B and crosses side 

rA of 7Z : so that the angle formed by BF and semi-line equals to ^ % ~H' A 1 Br , for 
each i = 1, . . . , m. These semi- lines split angle ABF into m—l angles, each of which 
is equal to 4 , ( -„^_ 1 - ) , since ABB = ir/4. Say </> = 4,(^-1) • Then, we place node uj at 
the intersection of semi-line U and rA, for each i = 1, ... ,771 (see Figllbj). In order 
to simplify the description of our approach, we denote by <z, the horizontal distance 
between two consecutive nodes u\ and u\ +1 , i = 1, . . . ,m — 1. 

We proceed by defining an additional bundle of m semi- lines, say £[,... ,£' m , that 
emanate from vertex A. More precisely, semi-line l[ emanates from vertex A and passes 
through the intersection of l m -i and rA (i.e., node u i m _ i ), for each i = 1, . . . , m (see 
Fig |lb[) . Let be the angle formed by two consecutive semi-lines V i and for each 
i = 1, . . . , m — 1. 

So far, we have managed to fix the position of the nodes of V\ only (along side 
FA of TV). Symmetrically, we define the position of the nodes of V2 along side AB of 
1Z. This only involves two additional bundles of semi-lines emanating from vertices 
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(a) A circular drawing of K, 



(b) A 2-layered drawing of K, 



Fig. 1: Illustrations of our constructions 



r and A. We now proceed to investigate some geometric properties of the proposed 
construction. 

Lemma 1. For each i = 1, 2, . . . , m — 1, it holds that a,_i < a^. 

Proof. By induction. For the base of the induction, we have to show that a± < a 2 . 
First observe that a\ = if tan and a\ + a-2 = H tan 20. Therefore: 



However, both tan0 and tan 20 are greater than zero, which immediately implies 
that a\ < a 2 . For the induction hypothesis, we assume that V k, k < m—1 it holds that 
<jfc_i < a k and we should prove that < <Xk+i- Obviously, a\ + . . . + a,k = H tan 
Based on Equation [1] and similarly to the base of the induction, we have: 

— ak+i = H tan <fi ■ (1 + tan (k + l)<fi ■ tan k(f>) 

— at = H tan cj> ■ (1 + tan (k — 1)0 ■ tan k<fi) 

In order to complete the proof, observe that (k + 1)0 > (k — 1)0. □ 

Lemma 2. For each i = 1, 2, . . . , m — 1, it holds that (j>' i _ 1 > 0J. 

Proof. By induction. For the base of the induction, we have to prove that 0' x > 0' 2 or 
equivalently that tan0^ > tan0 2 . It holds that tan0^ = a m _i/iJ and tan (0^ + 2 ) = 
(o m _i + a m ~2) / H . By combining these relationships with Equation Q] we have that 



which trivially holds due to Lemma [TJ For the induction hypothesis, we assume that 
V k, k < m — 1 it holds that <fi' k _ 1 > 4>' k and we have to show that (j>' k > 0' fc+1 - Observe 
that: 



a 2 = H (tan 20 — tan 0) = a± ■ (1 + tan 20 • tan 0) 




\ 3 2 
flm-2) + a m ~\ ' a rr 



a m -2 > 0, 
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- tan M — tan ^ 1 + ---+'^'t)~ tan ('^'i + ---+'^fc-i) _ Ha m k 

^fc ~~ l+tan(0i + ...+0i t )-tan(0i + ...+0' fc _ 1 ) ~~ if a + (a m _ 1 + ...+a m _ Jt )(a m _i + ...+a m _ fc+1 ) 

— tan 



fe+l^ ff 2 + (o m _i + ...+a m _ (fc+1) )(a m _i + ...+a m _ fc ) 

By Lemma [T]we have that (a m _i + . . . + a m _(fc_x)) > (a m -i + • • • + ^m-ffe+i)) an d 
■ff ' O'm—k > -ff ' o m -(Hi)' Therefore, tan</4 > tan^ , 1 . □ 

Lemma 3. Angle 4>' m _i is the smallest angle among all the angles formed in the 
drawing. 

Proof. From Lemma [51 it follows that angle 0j n _ 1 is the smallest angle among all 
<j>' i: i = 1, ... ,m — 1. Additionally, it is not difficult to see that angle 4>' m _\ is larger 
(but remains the smallest among all fy, i = l,...,m — 1), if the endpoint of the 
bundle (i.e., node u^) moves to any internal point (i.e., node uf) of side AB (see 
Fig llbjl . Therefore, 4>'m-i ^ s the smallest angle among all the angles formed by pairs 
of consecutive edges incident to any node of V^- Since m > n, the same holds for the 
nodes of V\. Therefore, <f>' m _i defines the angular resolution of the drawing. 

Consider now two crossing edges (refer to the bold, crossing dashed-edges of 
Fig llbj) . Their crossing defines (a) a pair of angles that are smaller than 90° and 
(b) another pair of angles that are larger than 90°. Obviously, only the acute angles 
participate in the computation of the crossing resolution (see angle 4> c in Fig |lb[) . 
However, in a complete bipartite graph the acute angles are always exterior to a tri- 
angle having two of its vertices on V\ and V2, respectively (refer to the gray-colored 
triangle of Fig llbj) . Therefore, the crossing resolution is always greater than the an- 
gular resolution, as desired. □ 

Lemma 4. It holds that 4>'m-i — \ - 

Proof. We equivalently prove that tan^_j > tan %. Using Equation^ we have that 
tan 2 < Iff ■ Therefore: 



tan (f> m _i > tan 



4> 



in 
H 



-1 a 2 + ...+a„ 



1 + 

1 T H H 

H 2 > (ai + . . . + Om-x)(02 

H-ai >0 



> 



2H 

dm-l) 



which obviously holds. □ 

Theorem 2. A complete bipartite graph K m>n admits a 2-layered drawing of total 
resolution 0( ma ^ m n} ). 

Proof. Immediately follows from Lemmata [3] and [4] □ 

Consider now the case where the nodes of the graph must have integer coordinates, 
i.e., we restrict ourselves on grid drawings. An interesting problem that arises in this 
case is the estimation of the total area occupied by the produced drawing. We will 
describe how we can modify the positions of the nodes produced by our algorithm in 
order to obey the grid constraints. Assume without loss of generality that L\ and £2 
are two horizontal lines, so that £2 coincides with y-axis and the drawing produced 



G 



by our algorithm has a\ = 1. Then, we can express the height of drawing r(K m ^ n ) as 
a function of <fi, as follows: 

ai = 1 <S=> tan <j> ■ H = 1 H = 1/ tan <f> 

Note that this drawing does not obey the grid constraints. To achieve this, we 
move the horizontal line C\ to the horizontal grid line immediately above it and each 
node of both V\ and V2 to the rightmost grid-point to its left. In this manner, we 
obtain a new drawing r'(K m n ), which is grid as desired. By Lemma [1] it follows 
that there are no two nodes sharing the same grid point, since a\ is slightly greater 
than one grid unit. Since neither horizontal line L\ nor any node of K m<n moves 
more than one unit of length, the total resolution of r'(K m , n ) is not asymptotically 
affected, and, in addition the height of the drawing is not significantly greater (i.e., 
asymptotically it remains the same). Based on the above, the area is bounded by 
cot 2 4> or equivalently by 1/ tan 2 <f>. By Equation 02 this is further bounded by l/<fi 2 - 
By Theorem[3J it holds that <fi = 0(1/ max{m, n}). Therefore, the total area occupied 
by the drawing is 0(max{m 2 , n 2 }). The following theorem summarizes this result. 

Theorem 3. A complete bipartite graph K m ^ n admits a 2-layered grid drawing of 
&{ r r) total resolution and 0(max{m 2 , n 2 }) area. 

4 A Force Directed Algorithm 

We present a force-directed algorithm that given a reasonably nice initial drawing, 
probably produced by a classical force-directed algorithm, results in a drawing of 
high total resolution. The algorithm reinforces the classical force-directed algorithm of 
Eades [S] with some additional forces exerted to the nodes of the graph. More precisely, 
these additional forces involve springs and some extra attractive or repulsive forces 
on nodes with degree greater than one and on end-nodes of edges that are involved 
in an edge crossing. This aims to ensure that the angles between incident edges and 
the angles formed by pairs of crossing edges will be as large as possible. The classical 
force-directed algorithm of Eades [5] models the nodes of the graph as electrically 
charged particles that repel each other, and its edges by springs in order to attract 
adjacent nodes. In our approach, we use only the attractive forces of the force-directed 
algorithm of Eades (denoted by Spring), which follow the formula: 

-^springO^Pt,) = ^spring ■ log — ■ p u pl, (u,v) £ E 

£ spring 

where C spr ing and Spring capture the stiffness and the natural length of the springs, 
respectively. Recall that p u p\ denotes the unit length vector from p u to p v . 

We first describe our approach for the case where two edges, say e = (it, v) and 
e' = (it', v'), are involved in a crossing. Let p c be their intersection point. W.l.o.g., we 
assume that it and u' are to the left of v and v', respectively, y u i < y u and y v < y v i , as 
in Fig(2l Let 9 vv i be the angle formed by the line segments p c p v and p c p v > in counter- 
clockwise order around it from p c p v to p c Pv' ■ 111 order to avoid confusion, we assume 
that 6 vv i = V ' V , i.e., we abuse the counter-clockwise measurement of the angles that 
would result in 6 vv i — 2ir — 8 V > V . Similarly, we define the remaining angles of FigJ^l 
Obviously, 9 VV > + 6 V > U = n. Ideally, we would like 9 vv i = 9 V > U = ^, i.e., e and e! form 
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Fig. 2: Forces applied on nodes in order to maximize the crossing resolution, (a) Springs on 
nodes involved in crossing, (b) Repelling or attractive forces based on the angles. 



a right angle crossing. As we will shortly see, the magnitude of the forces that we 
apply on the nodes u, u', v and v' depends on (a) the angles 9 VV > and 9 V > U and (b) the 
lengths of the line segments p c p u , p c p u > , p c Pv and p c p v ' ■ 

The physical model that describes our approach is illustrated in Figj^l Initially, 
for each pair of crossing edges at point p c , we place springs connecting consecutive 
nodes in the counter-clockwise order around p c , as in Fig |2al The magnitude of the 
forces due to these springs should capture our preference for right angles. Consider 
the spring connecting v and v'. The remaining ones are treated symmetrically. We 
set the natural length, say ^sp ring , of the spring connecting the nodes v and v' to be 

vHPc — P v \\ 2 + \ \P C ~ P'v\\ 2 - This quantity corresponds to the length of the line seg- 
ment that connects v and v' in the optimal case where 9 VV > = -|. So, in an equilibrium 
state of this model on a graph consisting only of e and e', e and e' will form a right 
angle. Concluding, the force on v due to the spring of v' is defined as follows: 



T CI0S (n n A - r CTOS lnr ^ Pv Pv ' 1 tT^T 5, 
•T spring \PviPv' ) — ^spring i0 B J^p PvPv' 



spring 

The remaining forces of Fig[2a] are denned similarly. Note that in the formula 
above, the constant C^°l ng is used to control the stiffness of the springs. 

Our preference for right angle crossings can be also captured using the angles 9 VV > 
and 9 V ' U (see Fig |2b[) . As in the previous case, we restrict our description on the angle 
formed by the line segments p c p v and p c p v > . Ideally, we would like to exert forces 
on the nodes v and v' such that: (i) when 9 VV > —> 0, the magnitude of the force is 
very large (in order to repel v and v'), and (ii) when 9 vv i — > the magnitude of 
the force is very small. A function, say / : M — > M., which captures this property is: 
f{9) = a . Having specified the magnitude of the forces, we set the direction of 
the force on v (due to v') to be perpendicular to the line that bisects the angle 9 VV > 
(refer to the dash-dotted line l vv < of Fig l2b| . or equivalently parallel to the unit length 
vector PeTp(Bsc(p c pl,p c pyi)). Recall that Perp and Bsc refer to the perpendicular and 
bisector vectors, respectively (see Section [2]). It is clear that if 9 VV > < the forces on 
v and v' should be repulsive (in order to enlarge the angle between them), otherwise 
attractive. This can be captured by the sign function. We conclude with the following 
formula which expresses the force on v due to v' . 

F£&b(Pv,Pv>) = C^°i e • sign{9 vv , - ^) • f{9 vv ,) ■ Perp(Bsc(p^,p c p„/)) 
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■Spring [Pu, , P, l+1 i U) ^anjilar [Pv, , > U) 

(a) (b) 

Fig. 3: Forces applied on nodes in order to maximize the angular resolution, (a) Springs on 
consecutive edges around u. (b) Repelling or attractive forces based on the angles. 



where constant C^°l e controls the strength of the force. Similarly, we define the 
remaining forces of Figl2bl 



Consider now a node u that is incident to d(u) edges, say erj = (u, vq), e\ = (it, v\). 
■ ■ ' ' e d(u)~i — ( u > v d(u)-i)- We assume that eo,ei, . . . ,&d(u)-i are consecutive in the 
counter-clockwise order around u in the drawing of the graph (see Fig l3a|) . Similarly 
to the case of two crossing edges, we proceed to connect the endpoints of consecutive 
edges around u by springs, as in Fig l3al In this case, the natural length of each spring, 
should capture our preference for angles equal to In order to achieve this, we 

proceed as follows: For each i = 0, 1, . . . , d(u) — 1, we set the natural length, say /s pring , 
of the spring connecting Vi with U(i + i)mod(ci(u))> to be: 

Spring = a/IMP + \\e( t+ i) m od(d(u))\\ 2 ~ 2 • ||e 4 | ■ \e [i+ i )mod{d ( u)) 1 1 • cos(2?r/d(w)) 

where ||e|| is used to denote the length of the edge e <E E in the drawing of the graph. 
The quantity ^ prlng corresponds to the length of the line segment that connects Vi with 
V(i+i) m od(d(u)) in the optimal case where the angle formed by a and e {i+1 ) mod ^ d{u)) is 
. Therefore, the spring forces between consecutive edges should follow the formula: 



T-angular/ „ . ql \ — ^angular , l|P»i_Pt >^+^ 

-'spring \PVi j Pv {i + 1 - }mod{d ^ u)) i ">j — "-"spring lu 6 p . 



-■PviPv 



(i + l)mod(d( U )) 



where the quantity C^pf^g 1 is a constant which captures the stiffness of the spring. 

Let now Qi be the angle formed by and e(i+i) mod (d(u)), measured in counter- 
clockwise direction from to e-(i+i) mo d(d(u)), i = 0, 1, ... , d(u) — 1. Similarly to the 
case of two crossing edges, we exert forces on V{ and i>(i+iw d(d(«)) perpendicular to 
the bisector of Qi, as illustrated in Fig|3bl However, in this case we need a magnitude 
function such that: (i) when Qi — > 0, the magnitude of the force is very large (in order 
to repel Vi and i>(j_|_i) TOO( j(d(u))); an d (ii) when Qi — > ^^y, the magnitude of the force 

is very small. Such a function, say j:Ix 7->K, is: g(Q; u) = — L. Having fully 
specified the forces applied on the endpoints of consecutive edges and their directions, 
we are now ready to provide the exact formulas that the forces follow: 
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F™ill* r {Pvi,Pv iz+1)mod{diu)) \u) = C^ e lar ■ sign(6i - • g(9i;u)- 

Perp(Bsc(p u p^ , PuPv ii+1)modWa) \ )) 

where C™|"g ar is a constant to control the strength of the force. In the work of Lin and 
Yen [17] , the above technique which applies large repelling forces perpendicular to the 
bisector of the angle, when the angle is small, is referred to as edge-edge repulsion. 
However, in their work, they use different metric to control the magnitude function. 
Of course, we could also use their metric, but we prefer the one reported above in 
order to maintain a uniform approach in both crossing and angular cases. In Section 
14. 1[ we provide an experimental comparison of these techniques. 

Note that by setting zero values to the constants C^°f ng , C^°* e or C^^ c , 

Wangle al \ our algorithm can be configured to maximize the angular, or the crossing 
resolution only, respectively. 

On each iteration, our algorithm computes three types of forces. Computing the 
attractive forces of the classical force-directed model among pairs of adjacent nodes 
of the graph requires O(E) time per iteration. The computation of the forces due to 
the edge crossings needs 0(E 2 ) time, assuming a straight forward algorithm that in 
0(E 2 ) time reports all pairwise crossing edges. Finally, the computation of the forces 
due to the angles between consecutive edges can be done in 0(E + Vd(G) \ogd(G)) 
time per iteration, where d{G) denotes the degree of the graph, since we first sort the 
incident edges of each node of the graph in cyclic order. Summarizing the above, each 
iteration of our algorithm takes 0{E 2 + Vd(G) \ogd(G)) time. 

The time complexity can be improved using standard techniques from compu- 
tational geometry [19j . If K is the number of pairwise-crossing edges, then the K 
intersections can be reported in 0(K + E log 2 Ej log log-E) time |191 pp.277], which 
leads to a total complexity 0(K + Elog 2 Ej log log E + Vd(G) logd(G)) per iteration. 



4.1 Experimental Results 

In this section, we present the results of the experimental evaluation of our algorithm. 
Apart from our algorithm, we have implemented the force directed algorithms of 
Eades [5] and Lin and Yen [17] • The implementations are in Java using the yFiles 
library (h ttp:/ /www.yworks.com| . The experiment was performed on a Linux machine 
with 2.00 GHz CPU and 2GB RAM using the Rome graphs (a collection of around 
11.500 graphs) obtained from graphdrawing.org. Fig[4] illustrates a drawing of a Rome 
graph with 99 nodes and 135 edges produced by our force directed algorithm. 

The experiment was performed as follows. First, each Rome graph was laid out us- 
ing the SmartOrganic layouter of yFiles. This layout was the input layout for all three 
algorithms, in order to speed up the experiment and overcome problems associated 
with local minimal traps especially in large graphs. If both the angular and the cross- 
ing resolution between two consecutive iterations of each algorithm were not improved 
more that 0.001 degrees, we assumed that the algorithm has converged and we did 
not proceed any more. The maximum number of iterations that an algorithm could 
perform in order to converge was 100.000. We note that the termination condition 
is quite strict and demands a large number of iterations. Our algorithm is evaluated 
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Fig. 4: A drawing of Rome graph grafol0129.99 consisting of 99 nodes and 135 edges with 
angular resolution 20.15° and crossing resolution 26.12°. 



as (a) Crossing-Only, (b) Angular-Only and (c) Mixed. The results are illustrated in 
Fig [5] and should be viewed in color. 

The Angular Resolution Maximization Problem: Refer to Fig |5al Our experimental 
analysis shows that our Angular-Only algorithm achieves, on average, better angular 
resolution. The angular resolution of our Mixed algorithm is almost equal, on average, 
to the one of Lin and Yen. Note that the algorithm of Lin and Yen, in contrast to ours, 
does not modify the embedding of the initial layout [T7] (i.e., it needs a close-to-final 
starting layout and improve on it). This explains why our Mixed algorithm achieves 
almost the same performance, in terms of angular resolution, as the one of Lin and 
Yen. In 59.76% of the graphs our Mixed algorithm yields a better solution compared 
to Lin- Yen's algorithm with an average improvement of 6.94°. 

The Crossing Resolution Maximization Problem: In Fig |5bl the data were filtered 
to depict only the results of non-planar drawings produced by the algorithms and 
avoid infinity values in the case of planar ones. It is clear that our Crossing-Only 
algorithm results in drawings with high crossing resolution. Our Crossing-Only al- 
gorithm performs better on large graphs compared to the algorithm of yFiles. The 
average improvement implied by our Crossing-Only algorithm is 13.63° w.r.t. the 
yFiles algorithm. 
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(a) Angular resolution results (b) Crossing resolution results 




(c) Total resolution results (d) Running time results 



Fig. 5: A visual presentation of our experimental results. The X-axis indicates the number 
of the nodes of the graph. In Fig.(a)-(c) the y-axis corresponds to the resolution 
measured in degrees, whereas in Fig.(d) to the running time measured in milliseconds. 



The Total Resolution Maximization Problem: Refer to Fig(Sc] This is the most 
important result of our experimental analysis. It indicates that our Mixed algorithm 
applied on graphs with more than 50 nodes constructs drawings with total resolution 
of 20 degrees, on average. Note that this is an achievement on average and under the 
particular termination condition discussed above, i.e., there is no guarantee, that it 
would be maintained indefinitely. An example is given in Fig0J 

Finally, Fig |5dl summarizes the running time performance of the algorithms. Our 
algorithm needs, on average, 7340 milliseconds and 1298 iterations to converge, whereas 
the one of Lin and Yen 8346 and 1952, respectively. Note that the time complexity 
of Lin- Yen's algorithm is better than ours. However, the termination condition takes 
into account the crossing improvement and therefore the algorithm of Lin and Yen 
needs more iterations to converge, which explains this contradiction. 



5 Conclusions 

In this paper, we introduced and studied the total resolution maximization problem. 
Of course, our work leaves several open problems. It would be interesting to try to 
identify other classes of graphs that admit optimal drawings. Even the case of planar 
graphs is of interest, as by allowing some edges to cross (say at large angles), we may 
improve the angular resolution and therefore the total resolution. 
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